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ABSTRACT 

The Extreme Ultraviolet Explorer (EUVE) satellite accumulated ~ one million seconds of 
public data between 1994 and 1997 for the BL Lacertae object Markarian 421. This is the second 
of two papers in which we present the results of spectral and temporal analysis of this EUVE 
data set. We analyze in the present paper the imaging data by means of power spectrum and 
structure function techniques, while the spectral analysis is presented in a companion paper. We 
find for Mrk 421 a power spectrum with slope —2.14 ± 0.28 with a break at ~ 3 days. This is 
the first time that a break in the power spectrum of a BL Lacertae object has been found. We 
also find evidence of non-stationarity for Mrk 421 EUV emission. 

Subject headings: galaxies: nuclei - galaxies: active - galaxies: individual (Mrk421) 



1. Introduction 

BL Lacertae objects are an extreme subclass 
of Active Galactic Nuclei (AGNs) emitting highly 
variable non-thermal radiation over 20 decades 
of frequency from radio to TeV energies. The 
mechanism responsible for the production of ra- 
diation over such a wide range is thought to be 
synchrotron emission followed by inverse Comp- 
ton scattering from energetic electrons. Some of 
the extreme properties of BL Lac objects, such 
as rapid variability, high gamma-ray luminosities 
and superluminal motion (e.g. Urry & Padovani 
1995), are explained invoking relativistic beaming 
of the emitted radiation along the jet axis which 
is taken to be oriented at small angles to the line 
of sight. 

Mrk 421 is one of the closest (z=0.0308) and 
best studied BL Lacertae objects. It is one of the 
strongest TeV sources and one of the brightest ex- 
tragalactic EUV sources. The EUV emission of 
Mrk 421 during quiescence is believed to be the 



peak of a synchrotron spectrum. Mrk 421 was 
monitored by the EUVE satellite for four years 
(1994,1995,1996 and 1998) as part of multiwave- 
length campaigns and it was serendipitously ob- 
served in 1995 and in 1997 with the EUVE pho- 
tometers for a total of one million seconds, the 
best sampled light curve for any BL Lac object at 
high energies 1 . Table 1 summarizes all the EUVE 
observations. 

The main purpose of the present paper is to an- 
alyze in an homogeneous way the 1100 ks of EUVE 
public imaging data for Mrk 421 (less than half of 
it has been analyzed to date, Fruscione et al. 1996, 
Kartje et al. 1997) in order to perform the first de- 
tailed variability analysis (both in terms of power 
spectrum and structure function analysis) for a BL 
Lacertae object in the EUV. In general, variabil- 

As in paper I, wc are not considering the most recent ob- 
servation (April 1998) since the data are presently being 
analyzed by the original proposers (Marshall et al. in prepa- 
ration) 



ity studies by means of power spectrum analysis 
can reveal characteristic time scales of the emis- 
sion mechanism (e.g. periodicities) and can con- 
strain physical models. Our analysis reveals the 
existence of a characteristic timescale for Mrk 421 
that can constrain current models. 

The paper is organized as follows: all the EUVE 
lightcurves of Mrk 421 collected from 1994 to 1997 
are presented in Section 2. In Sections 3 and 4 we 
present the results from the power spectrum and 
structure function analysis, respectively. In Sec- 
tion 5 we discuss those results. Finally in Section 
6 we give some conclusionary remarks. For the 
energy spectrum data and analysis we refer the 
reader to the companion paper Cagnoni & Frus- 
cione (2000, hereafter paper i). 

2. Observations 

A general description of the Extreme Ultravio- 
let Explorer satellite (EUVE) is given in paper I. 
The photometric instruments on board EUVE are 
the Deep Survey/Spectrometer, (DS/S, covering 
the 68-178 A band for photometry; Welsh et al. 
1990) used for the deep EUV survey along the 
ecliptic and for pointed spectroscopic and imag- 
ing observations, and three co-aligned photome- 
ters (the "scanners"), observing in four different 
bandpasses (from 58 to 740 A) and mounted or- 
thogonally to the DS/S. The photometers were 
used to carry out the EUVE all-sky survey and 
are now employed in the "Right Angle Program" 
(i.e. the observation of serendipitous targets while 
the DS/S is conducting the primary science obser- 
vation). Figure 1 shows a comparison of the DS/S 
and scanner effective areas. 

Mrk 421 was observed several times from 1994 
to 1997: four times with the DS/S and twice with 
the scanners; Table 1 summarizes these observa- 
tions. Because of interstellar medium absorption 
along the line of sight (Nh = 1.45 x 10 20 cm~ 2 , 
Elvis, Wilkes & Lockman 1989), Mrk 421 was de- 
tected only in the Lexan/B filter (the shortest 
wavelength bandpass) both in the scanners (58- 
174 A) and in the DS/S (67-178 A). 

2.1. Serendipitous sources 

Two faint serendipitous sources are clearly vis- 
ible next to Mrk 421 in the DS image. We 
identified the closest one (~ 2' from Mrk 421 



) with 2 unresolved stars in a double system: 
HD 95934 (M v = 6.00, Spectral type A3III-IV) 
and BD+39 2414B(My = 12.6). The other one at 
<~ 4' from Mrk 421 was identified with another star 
in double system: HD 95976 (M v = 7.4, Spectral 
type FII). 

The DS 50% energy radius is about 0.3' up to 
0.3° off-axis. In scanner A, this half energy radius 
is <~ 3' at <~ 1.5 — 2° off-axis (where the sources lied 
in the 1995 and 1997 observations) and the sources 
are not resolved. During the analysis of DS data 
we excluded a circle centered on HD 95934; this 
excluded only a negligible fraction of counts from 
Mrk 421 . 

2.2. Light curves 

We extracted the light curves from the DS or 
Scanner time-ordered event lists using the EUVE 
Guest Observer Center software (IRAF/EUV 
package) and other IRAF timing tasks adapted 
for EUVE data. All the lightcurves were binned 
over one average EUVE orbit (5544 s), but effec- 
tive exposure times per bin vary due to several 
correction factors different case by case (e.g. in- 
strumental deadtime, telemetry saturation, high 
particle background, dead spot). 

Unless described otherwise, we chose the ex- 
traction region for the source as a circle with a 
2.5' radius excluding a circle with a 0.67' radius 
centered on the serendipitous source. The back- 
ground region is an annulus with inner and outer 
radii equal to ~ 4' and 8' respectively which avoids 
the other serendipitous source. 

In order to compare scanner and DS data, we 
had to normalize the count rates. Count rates de- 
pend on the effective area of the telescope and on 
the spectral shape of the source. We convolved a 
sample of input spectra (absorbed power law with 
slopes as in Table 2 of paper i) with the DS and 
scanner effective areas and defined the normaliza- 
tion factor as the ratio of peak of the derived count 
rates (see Figure 2 as an example for a = 2). The 
average normalization factor is about 1.4±0.3 (de- 
pending on the input spectrum). We used this 
value throughout the paper and verified that the 
results in the power spectrum analysis are all con- 
sistent within the errors of this value. 

The combined light curve of all the observations 
is presented in Figure 3. In this section we will 
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give a qualitative description of each lightcurve 
in terms of flares and sub-flares: a general ris- 
ing/decreasing trend will be addressed as a "flare" , 
while faster variations superimposed on the gen- 
eral trend will be referred to as sub- flares. This 
description is useful in understanding the results 
from the power spectrum and structure function 
analysis. 

2.2.1. 1994 DS light curve 

Figure 4 shows the 1994 light curve, which 
was obtained eliminating time intervals when the 
source crossed the dead spot 2 (about 1440 s). 

The source dimmed for about 3 days, bright- 
ened again for another ~ 3.5 days, until it reached 
a level of ~ 0.2 counts/sec and stayed there for 
the remaining of the observation (another ~ 3.5 
days). This light curve behavior could be sam- 
pling the end of a flare and the beginning of a new 
one. Interestingly, when the source increases its 
flux, after MJD 49447, (i.e. possibly during the 
beginning of a new flare) smaller "sub-flares" are 
superimposed on it; this trend is also present in the 
1995 scanner lightcurve (see 2.2.2 and Figure 5). 

2.2.2. 1995 scanner A light curve 

During this scanner exposure optimized for the 
primary scientific DS/S observation of HD 41511, 
Mrk 421 partially falls in the region usually dis- 
carded by the standard analysis pipeline because 
of the obscuration due to the filter support and 
because of the high particle background close to 
the detector edge. We extracted the lightcurve 
only from the "clean" region of the detector and 
we corrected the count rate for the fraction of 
the count lost (<~10%). The background region 
was chosen in the same detector quadrant as the 
source and large enough to match the ratio back- 
ground/source (~ 7.8) region used for the other 
lightcurves. The lightcurve is plotted in Figure 5; 
the open circles represent the count rate mea- 
sured in scanner A, while the filled circles are the 
count rate normalized to the DS; this lightcurve 
has as behavior similar to the 1994 DS lightcurve: 

2 A small region of reduced gain and detector quantum ef- 
ficiency near the center of the DS detector, caused by the 
observation of the very bright EUV source HZ 43 (Sirk 
1994) 



when the source brightens, after MJD 49754, "sub- 
flares" appear above the main trend. 

2.2.3. 1995 DS light curve 

During this DS/S observation, the source was 
pointed 0.3° off axis both to avoid the dead spot 
and to extend the spectrum toward shorter wave- 
lengths (see paper i). The light curve is shown in 
Figure 6. 

The behavior of the source during the DS 1995 
observation is complex: at the beginning Mrk 421 
brightens by a factor of ~ 65% in less than 2 days 
and this flare, which lasts for about 3 days, is well 
correlated to a TeV/X-Ray flare (Buckley et al. 
1996). Right after the first one, there is a sec- 
ond flare lasting for ~ 6 days which can also be 
seen as the superposition of two separate flares 
(the second one starting at ~ MJD 49838). After 
MJD 49845.5, the source rises from ~ 0.20 counts 
s _1 , up to <~ 0.28 counts s _1 , and then decays to 
<~ 0.23 counts s _1 , within ~ 3 days. The overall 
variability of the source is about a factor 3. 

2.24. 1996 DS light curves 

In 1996 Mrk 421 was observed few days before 
and few days after the detection of a large TeV 
flare on May 7th 1996 (Zweerink et al. 1997). The 
lightcurves are shown in Figure 7. This lightcurve 
looks different from all the other ones because no 
slow variations seem to be present; the source ap- 
pears to be at a constant level with lots of small 
amplitude, fast variations superimposed on it; i.e., 
the light curve appears to be the sum of flares 
which live a short time, and rise or decay faster 
when compared with the flares seen in the other 
light curves (see for example the two flare like 
events after MJD=50200). 

2.2.5. 1991 scanner A light curve 

In 1997 Mrk 421 was observed for 4 days with 
the scanner A telescope during the DS pointed 
observation of T Tauri; the source was detected 
within the useful area of the scanner A field of 
view allowing an easier data reduction than the 
previous scanner observation. The light curve is 
shown in Figure 8: the open circles represent the 
scanner data while the filled circles correspond to 
the count rate normalized to the DS. This is the 
shortest lightcurve and its shape suggests that we 
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may be looking at the tail of a flare and at the 
beginning of a new one. 

3. Power Spectrum Analysis 

In order to quantify the variability seen in the 
EUVE lightcurves of Mrk 421 , we computed their 
pcriodogram as: 



*M = *]e- a * VtU Y C, (1) 

where: x, At and N are the mean value, bin 
size and number points of each light curve respec- 
tively and Vi = i/(NAt), i = 1, 2, (N/2) - 1 
(e.g. Deeming 1975). The constant C in equa- 
tion (1) represents the power level due to Poisson 
statistics, C = Ata 2 nmse = (At/N)J2? =1 err(i) 2 
(err(i) is the error of each light curve point). The 
subtraction of this constant level is necessary since 
we are interested in the source power spectrum 
only and the EUVE light curves were obtained 
with two different instruments, hence the Pois- 
son statistics will be different in them. The pe- 
riodogram calculated in this way has the units 
of (counts/sec) 2 /Hz, and its integral over positive 
frequencies is equal to half the light curve variance. 

It is common to average the periodogram es- 
timates in order to reduce the variance and ap- 
proximate their probability distribution function 
to a Gaussian. Furthermore, it is better to esti- 
mate the logarithm of the power spectrum instead 
of the power spectrum itself, since the logarithmic 
transformation brings the distribution function of 
the power spectrum estimates closer to a Gaussian 
(Papadakis & Lawrence, 1993). Therefore, we cal- 
culated the logarithm of the binned periodogram 
estimates, 

l O g 10 [P(^)]=log 10 [^A^], (2) 

(where i/j = J^j Vj/m, and i/j = j/(NAt)) and 
used it as our estimate of the logarithm of the 
power spectrum. The index j varies over m con- 
secutive periodogram estimates and the index i 
varies over the [(N/2) — l]/m subsets into which we 
have grouped the periodogram estimates. Taking 
into account the propagation of errors (e.g. Bev- 



ington 1969), the error on the logarithm of the 
power spectrum is 

<Klogu#(^)]) = logio( e ) cr p(„ l) / i3 (^)' ( 3 ) 

where is the square root of the variance 

of the I(vjy s around the mean value, P(vi), at 
each bin. 

3.1. Power spectra of the orbit-binned 
lightcurves 

We estimated the power spectrum of the DS 
and scanner A light curves binned over one av- 
erage EUVE orbit (5544 s) shown in Figures 3 - 
8. To avoid the gap in 1995 DS observation (Fig- 
ure 6) we divided it into 2 parts and performed a 
power spectrum analysis on each one of them. Us- 
ing equations 1 and 2, we estimated the logarithm 
of the power spectrum with m = 10. A larger 
value for m would bring the distribution function 
of the power spectrum estimates closer to a Gaus- 
sian (Papadakis & Lawrence 1993). However this 
small value of m was necessary in order to have 
a reasonable frequency resolution, i.e. in order 
to have the maximum possible number of power 
spectral points above the Poisson noise level. 

The power spectra of all the light curves, nor- 
malized to the square mean value of the respective 
light curve, are shown in Figure 9. The normalisa- 
tion is necessary if one wants to compare the power 
spectrum amplitudes, since the light curves were 
obtained from two different instruments. All the 
power spectra show a similar red noise behavior at 
frequencies < 10~ 4 - 3 Hz (i.e. all of them rise loga- 
rithmically towards lower frequencies in a similar 
way), except for DS 1996 observation which seems 
to flatten for v < 1CT 5 Hz. 

We fitted the DS power spectra with a sim- 
ple power law model of the form P(v) = A x 
(1CP 5 ^)~ Q (the scanner A power spectra had too 
few points to fit a model to them). In this model, 
the normalization A is equal to the value of P(v) 
at v = 10~ 5 Hz. The model fit was done using 
the power spectra before normalizing to the mean 
square value. The best model fit results are listed 
in Table 2. The model gives a good fit to all the 
DS power spectra, but, due to the small number 
of points, the errors on the model parameters are 
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large. The best fit values of A are within the er- 
rors for all the DS spectra. However, the best fit 
slope values show a larger scatter. In particular, 
the difference between the DS 1995 (part 1) and 
the DS 1996 best fit slope values is 2.7 ±1.0. This 
~ 3cr difference indicates that the shape of the 
power spectrum docs not remain constant with 
time, i.e. that the EUV emission of Mrk 421 is 
not stationary 3 . However, looking at DS 1996 
power spectrum (Figure 9), there is a possibility 
that the flatter slope obtained from the fit is due 
to a flattening around v = 10~ 5 Hz. If this is the 
case, a broken power law model should give a bet- 
ter fit to the data. However, with just 5 points, 
it is not possible to investigate if a broken power 
improves the fit to the data in a statistically sig- 
nificant way (e.g. by the use of an F-test). On the 
other hand, perhaps the flattening of the power 
spectrum is a spurious effect introduced by the 
small binning used (m = 10). In an attempt of 
ruling out this possibility we used a larger value 
(m — 20). FigurelO compares the results: open 
circles for m = 10 (as in Figure 9) and filled cir- 
cles for m = 20. Both power spectra look similar, 
and the lack of frequency resolution in the latter 
case cannot allow us to reach a firm conclusion. 

Despite this possible non-stationarity, it is use- 
ful to combine the periodograms of all the DS and 
scanner A light curves in order to increase the 
number of power spectrum estimates (and hence 
the frequency resolution) . We first normalized the 
data sets to the square mean value of the respec- 
tive light curve and then combined them into one 
dataset by sorting them in order of increasing fre- 
quency. Then we grouped this "combined" peri- 
odogram into bins of size m = 20 since the larger 
number of estimates allows us to use a larger bin 
size value. Figure 11 shows a plot of the result- 
ing power spectrum (filled circles). Note that the 
combined power spectrum has the units of Hz -1 . 
This spectrum can be considered as the "average" 
EUVE power spectrum of the source. A power 
law model gives a good fit to this average power 
spectrum (solid line in Figure 11; Table 2) with a 
best fit slope value of 2.14 ± 0.28. 

A lightcurve is said to be stationary if its mean, variance 
and power spectrum remain constant in time (e.g. Chat- 
field 1989); sec Section 5.4 for further details 



3.2. Power spectra of the 500 s binned 
lightcurves 

Due to the indication of a possible break in the 
power spectrum and/or a possible non-stationarity 
in the EUVE light curve of MRK 421, we decided 
to investigate further the power spectra of the DS 
1995 and 1996 light curves. Our main aim was to 
determine the high frequency shape of the power 
spectrum as well as possible. For this reason, we 
used a ~ 10 times finer binning, with a bin size 
of 500 s, for the light curves; we were thus able to 
extend the power spectrum to higher frequencies 
and to use a larger value for m (m = 20 in equation 
2) to bring the distribution function of the power 
spectrum estimates closer to a Gaussian. 

Figure 12 shows the DS 1995 and 1996 power 
spectra estimated with equation (1), without sub- 
tracting the expected Poisson noise level, and 
binned as in equation (2) with m = 20; the long 
dashed lines in these plots indicate the Poisson 
noise power level C. Both power spectra show 
a large peak at ~ 1.8 x 10~ 4 Hz corresponding 
to a period of 5556 s, i.e. the EUVE orbital pe- 
riod. The next strongest peak is at a harmonic of 
that frequency. These peaks are due to the uneven 
sampling of the 500 s binned light curves. 

Apart from those peaks, the two power spectra 
appear to have a similar shape: they rise towards 
low frequencies, and they are flat at higher fre- 
quencies. So, it seems as if the use of a smaller 
bin size in the light curves has not helped us much. 
The power spectra, at most frequencies, appear to 
be flat and dominated by the Poisson noise power. 
However, while the DS 1995 estimates at high fre- 
quencies are close to the expected Poisson power 
level, the DS 1996 power spectrum is higher than 
this level. 

We fitted both spectra with a modified power 
law model with one of the slopes fixed at 2.14 (see 
below). Due to the effects of the window func- 
tion, it is not possible to fit the logarithm of the 
model directly to the power spectrum calculated 
from the 500 sec binned light curves. Instead, we 
convolved the model, let us say M(u), with the 
window function at each frequency fj, 

W"0 = -jj / M{v)W{ Vl - v)dv. (4) 

" J — oo 
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The function Imdiiyi) takes account of the win- 
dow function effects. We call I m di{vi) the "con- 
volved model" . 

As M(v), we used the following function, 

MM = P^W + C ' (5) 

which has a power law form with a slope of 2a 
at frequencies >> fb/2n, and flattens to zero slope 
at frequencies << fb/2ir (as before, C represents 
the power level due to the Poisson counting statis- 
tics). Note that as M(v) we did not use a simple 
power law model because in this case it is not easy 
to estimate the convolution model (equation 4) an- 
alytically (a simple power law model cannot be de- 
fined at v = 0) and because there might indeed be 
a break, as suggested by DS 1996 power spectrum 
(Figure 9). 

Then, we averaged the convolved model I m di (f p ) 
over groups of size m = 20 exactly as we did when 
estimating the power spectrum. The logarithm of 
the convolved model power spectrum becomes: 

logiot-PmcM^i)] = logio fi — — ]■ (6) 

m 

We can now fit log 10 [.P m( «(i/j)] to the estimated 
power spectrum, log 10 [P(fj)]. In this way, we are 
taking into account both the window function and 
any effect that the binning procedure may have on 
the power spectrum estimation. 

For the fit we used standard % 2 statistics, we 
kept 2a fixed to the value 2.14, and we did not 
allow C to vary but rather used the value found 
as described at the beginning of Section 3. Table 3 
summarizes the best fit results for each observa- 
tion; the errors quoted represent the 90% confi- 
dence levels assuming that A and fb are correlated. 
By far the largest x 2 value is for the DS 1996 
model fit. However, it is not straightforward to 
judge the goodness of fit using the these % 2 values. 
Due to the uneven sampling of the light curves, 
the periodogram estimates are no longer indepen- 
dent. Therefore, the errors on log 10 [-P(^i)] ma y 
be larger than those estimated (and hence the x 2 
values smaller). At the same time, depending on 
the amount of correlation between the /(i^)'s, the 
log 10 [P(fj)] estimates may be correlated as well. 
As a result, the number of degrees of freedom may 
be smaller than those listed in Table 3. 



The best model fits are also plotted in Figure 
12 with a solid line; these plots show that the 
model agrees well with the DS 1995 power spec- 
trum, but not with the DS 1996 (as expected from 
the large x 2 value). We fitted again the DS 1996 
power spectrum with the model defined in equa- 
tion (5) but we let the slope free parameter 
this time. The best fit results are listed in Table 3 
(values in a parenthesis), and the best fit model 
is also plotted in Figure 12 (b) (dotted line). As 
this figure shows, the new model fit agrees bet- 
ter with the power spectrum. The best fit slope is 
now 1.1±0.1 significantly smaller than the average 
slope of 2.14 ± 0.28. The x 2 value is also reduced 
by 74.9. The ratio of the reduced x 2 values for the 
two fits is 2.1 (note that this ratio does not depend 
on whether the errors on the power spectrum are 
underestimated or not). Using the number of de- 
grees of freedom listed in Table 3 and the F-test, 
we find that the improvement on the fit when we 
let the slope to be a free variable is highly signif- 
icant (the probability of a significantly better fit 
when a = 1.1 is larger than 99.5%). The probabil- 
ity remains larger than 95%, even if the degree of 
freedoms listed in Table 3 are reduced by a factor 
of two. 

Furthermore, to check for any other possible 
indication of non-stationarity, we performed a fit 
to the power spectra of all the 500 sec binned light 
curves (including the scanner A light curves) using 
the model defined in equation (5) with 2a = 2.14 
. The best fit model results are listed in Table 3, 
and the best fit model, together with the power 
spectra, arc plotted in Figure 13 (the dotted lines 
in this figure show the Poisson noise power level). 
This figure shows that the model gives a good fit 
to all the power spectra. 

We conclude that all the EUVE power spectra 
of Mrk 421 are consistent with a power law model 
with a slope of <~ 2.1. The only exception is the 
power spectrum of the DS 1996, which shows a 
significantly flatter slope of 1.1 ± 0.1. The clear 
change in slope in the 1996 data shows that the 
EUV emission of Mrk 421 is not statistically sta- 
tionary (see Section 5.4). 

As hinted by Figure 9 for DS 1996 (open cir- 
cles), from the fit of the 500 s binned power spec- 
tra we found indication for the presence of a low 
frequency turn over in the EUVE power spectrum 
of Mrk 421 below which the power spectrum flat- 
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tens. This finding, if confirmed, would be the first 
detection of a break in the power spectrum of a BL 
Lac object. Since we cannot estimate the power 
spectrum at frequencies lower than ~ 10~ 5 Hz, 
using the value m — 20, we cannot be certain of 
the slope below the turn over frequency. All we 
can say with the present data is that the power 
spectrum should show a change in its slope at /;,. 
Furthermore the value of this knee frequency is 
not constrained very well from our data. Using 
the best fit results from the DS spectra only (as 
they have the largest number of data points), we 
get a straight mean value of /& = 4 x 10 ~ 6 Hz. By 
simply combining the overlapping 90% confidence 
limits on the best values of we find that with 
probability larger than 68%, fb should be between 
7 x 10~ 7 Hz and 8 x 10~ 6 Hz. 

From a physical point of view, it is natural to 
expect that low frequency knees should be present 
in the power spectra of AGN. Lack of this turnover 
would imply a physical system with infinite vari- 
ance, a system whose properties change contin- 
uously with time. The frequency turnover is not 
evident in the combined EUVE power spectrum of 
Mrk 421 (Figure 11), because that spectrum does 
not extend to frequencies below ~ 10 -5 5 Hz. Its 
presence is implied only when we take proper ac- 
count of the window function effects in the model 
fitting of the 500 s light curve power spectra. If 
the frequency turn over were located at a much 
lower frequency, more power would be transfered 
to higher frequencies. The resulting power spec- 
trum would have a higher normalization than the 
observed one. To put in other words, the ampli- 
tude of the variations in the light curves would 
be larger, i.e. they should have a higher variance. 
For example let us assume that the source power 
spectrum is well describe by equation (5) and that 
the break is at 10~ 7 Hz. Using equations (4) and 
(6) we can predict how the power spectrum of the 
DS 1995 light curve should be in this case (dot- 
ted line in Figure 12a). As this Figure shows, the 
power spectrum should have a much larger ampli- 
tude than the one that is actually observed. 

Due to the significance of such a feature in the 
power spectrum of the source, we searched further 
evidence of its existence by means of a structure 
analysis of the light curves. We present below the 
results from this analysis of the EUVE light curves 
which confirm the existence of this turn over in the 



power spectrum of Mrk 421 . 

4. Structure Function Analysis 

The Structure Function (SF hereafter) is de- 
fined as, 

SF(r) = 2a 2 -R(r), (7) 

where a and R{t) are the variance and auto- 
covariance at lag r of the light curve (equation A4 
in Simonetti, Cordes & Heeschen, 1985). For a 
given light curve, it is usually estimated by, 

SF(r)^^ N ^[x( l + r)- X (i)}^ (8) 

(equation A8 in Simonetti, Cordes & Heeschen, 
1985). The SF gives a measure of the mean dif- 
ference in the flux of two light curve points, as a 
function of their time separation (lag) r. There is 
a simple correspondence between the SF and the 
power spectrum: when P{v) cx v~ a then SF(t) cx 
t" -1 . Furthermore, if the power spectrum flattens 
below a frequency fb, then SF(t) ~ 2c 2 in the 
limit T»l/f b . 

The main advantage of the structure function 
analysis is that the SF estimation is not affected 
by gaps of missing points in the light curve (it is 
easy to show using equation (8), that < SF(t) >= 
SF(t), irrespective of the light curve sampling 
pattern). However, it is difficult to estimate 
the errors on SF(t) because, as with estimates 
of other time domain quantities (eg auto, cross- 
correlations), the SF estimates at different lags 
are heavily correlated (much more so than the pe- 
riodogram points at different frequencies, even in 
the case of unevenly sampled data). 

To estimate the SF of MRK 421, we used to- 
gether all the orbit-binned DS and scanner A light 
curves and constructed a single, long EUVE light 
curve of the source (in order to use the scanner 
data into an overall light curve, we normalized 
them to the DS expected count rate as explained 
in Section 2.2). For each pair of data points in the 
light curve, x(tj) and x(tj), we computed their 
time difference, Atij = t» — tj, and their flux dif- 
ference, FDij = x(ti) — x(tj). Let us assume 
that we have chosen a lag bin of size At, and 
that there are n points with Atij in the interval 
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[r — At/2, t + At/2]. Wc used the average value 
of the n FDij measurements, ^FDij/n, as an 
estimate of the structure function, SF(t'), at the 
lag t' = J2 At i,j/n. 

For the overall EUVE light curve (Figure 3), 
we estimated the SF twice. The first time we 
used At = 1 day, and the second time At = 1 
month. The first bin value was used to estimate 
the SF from lag ~ 1 day up to lag ~ 18 days. 
We cannot estimate the SF at higher lags with 
At = 1 day since the individuali?[/V7? light curves 
of Mrk 421 have a maximum length of ~ 20 days. 
Consequently, we used the larger bin size (At = 1 
month) to estimate the SF at lags > 2 months and 
up to ~ 2.5 years (ie the time interval between the 
first and last EUVE light curve of Mrk 421 ). 

The combined short and long bin estimated SF 
is plotted in Figure 14 (filled circles for At = 1 
day and filled squares for At = 1 month). The 
estimated SF rises roughly logarithmically with r 
at lags < 20 days. However, at lags > 60 days it 
is relatively flat. This flattening cannot be due 
to the uneven sampling of the light curve, but 
it is exactly what we would expect if there is a 
low frequency turn over in the power spectrum of 
Mrk 421 . 

We examine below whether the estimated SF 
shown in Figure 14 is consistent with the assump- 
tion that the true power spectrum of Mrk 421 is 
similar to the model defined by equation (5) with 
2a = 2.14 and f b = 4 x 10~ 6 Hz. 

To answer this question, we performed a nu- 
merical experiment: we created one set of 1000 
simulated light curves using the method described 
in appendix C of Papadakis & Lawrence (1995). 
Each simulated light curve has a power spectrum 
with slope of 2.14, a knee frequency at 4x 10~ 6 Hz, 
and variance as the EUVE light curve. For each 
one of the simulated light curves, we estimated 
the SF exactly in the same way as we did with 
the EUVE light curve. Then we computed the ex- 
pected mean value, < SF sim (T) >, and variance, 
a <SF ■ >( T )> °f au the simulated SFs at each r. 
The results are plotted in Figure 14 (open cir- 
cles and open squares for the mean short and long 
bin simulated SF). The "errors" on the points are 
equal to ±cr < sF 3im >('r). Note that these errors 
would represent the 68% of the process realiza- 
tions in the case of a Gaussian distribution, but, 



since the actual SF are not normally distributed, 
they can only give us an idea of the real la errors 
of the actual SF and cannot be used in the normal 
X estimation. In order to quantify the compar- 
ison between log[< SF sim (T) >] and log[5F(T)] 
in a statistical way, we estimated the sum of their 
squared differences, i.e. X^0°g[< SF sim (Ti) > 
} — \og[SF(ri)]} 2 , for all the for which we have 
estimated the SF. We refer to this sum as "x 2 '\ 
although its definition is different from the tra- 
ditional definition of x 2 ■ 111 our case, x 2 = 4. 
Based on the distribution of the x 2 values for the 
1000 simulated light curves, the probability that 
the EUVE light curve of Mrk 421 is a realization 
of a process with a power spectrum with a knee 
frequency at 4 x 10~ 6 Hz is ~ 40 %. 

Finally to check whether Mrk 421 structure 
function is consistent with a break at time scales 
longer than the range suggested by the fit of the 
power spectra, we performed a second numerical 
experiment: we created a new set of 1000 sim- 
ulated light curves, with 2a = 2.14, variance as 
the EUVE light curve variance, and a knee fre- 
quency at 10 Hz. As above, for each simulated 
light curve we estimated the SF, and then com- 
puted < SF S i rn (T) >, and <J < sF aim > (t). The re- 
sults for the mean short and long bin simulated 
SF are plotted in Figure 14 (dotted triangles). At 
large lag values the mean simulated SF agrees well 
with the SF of Mrk 421 . This is not surprising 
since, as we mentioned at the beginning of this 
Section, we expect that SF ~ 2a 1 in the limit 
when t is larger than l/fb (10 7 s in this case). 
However, the agreement between the mean simu- 
lated SF and the Mrk 421 SF at short lags is not 
good. The x 2 this time is equal to 183. Based on 
the distribution of the x 2 values for the 1000 simu- 
lated light curves, the probability that the EUVE 
light curve of the Mrk 421 is a realization of a pro- 
cess with a power spectrum with a knee frequency 
at 10~ 7 Hz is < 0.1 %. 

We conclude that the SF of Mrk 421 (Figure 14) 
is entirely consistent with the results from the 
power spectrum analysis (Section 3); the slope of 
the Mrk 421 power spectrum is <~ —2.1 and there 
is a knee frequency at <~ 4 x 10 -6 Hz, below which 
the spectrum flattens; in fact we found that it is 
highly unlikely that this knee frequency will be 
smaller than 10 -7 Hz. 
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5. Discussion 

5.1. Summary of results 

All the EUVE light curves of Mrk 421 show sig- 
nificant variations (Figures 4-8). One interpreta- 
tion is that they are the result of the superposition 
of flares whose duration, roughly estimated from 
the same figures, is about 3 days; as an example 
we quote the 1994 lightcurve (Figure 4) where the 
source flux declines for about 3 days, rises for an- 
other ~ 3.5 days and stays steady to the end of 
the observation (~ 3.5 days). We investigated the 
variability present in the light curves by means of 
a power spectrum analysis and our results are as 
follows: 

• on average, the EUVE power spectrum of 
Mrk 421 has a power law form with a slope 
-2.14 ±0.28; 

• there is a low frequency turn over at 4 x 1CT 6 
Hz (i.e. ~ 3 days, the 68% confidence range 
is 7 x 1CT 7 Hz - 8 x 1CT 6 Hz). We cannot 
determine the shape of the power spectrum 
below this turn over, nevertheless the data 
so far is consistent with a flattening of the 
power spectrum below this frequency to a 
slope of zero; 

• there is a strong indication that the power 
spectrum does not remain the same, i.e. that 
the EUVE emission of Mrk 421 is not sta- 
tionary. The DS 1996 observation shows 
a power spectrum with a slope —1.1 ± 0.2, 
which is much flatter than the average slope. 
However, this observation is again consistent 
with a flattening below the same frequency 
of 4 x 10~ 6 Hz. 

Furthermore, we performed a structure func- 
tion analysis of the total EUVE light curve. The 
structure function of the source flattens at time 
scales > 10 6 s. Based on numerical simulations, 
we found that the Mrk 421 SF is entirely consis- 
tent with a 2.14 power spectrum which flattens to 
a slope of zero below 4 x 10~ 6 Hz. 

We will discuss each of these results in more 
detail in the following subsections. 



5.2. The power spectrum slope 

Our analysis shows that, on average, the EUVE 
power spectrum of Mrk 421 (up to ~ 10 -4 5 Hz) 
has a power law form with a slope of —2.14 ±0.28. 
No periodicities or quasi-periodicities could be 
found. To the best of our knowledge, this is the 
first time that the power spectrum of a radio-loud 
AGN has been estimated in the EUV. In fact, even 
in the X-rays, the power spectra of radio-loud ob- 
jects are not well known. The only exception is 
PKS 2155-304. Its X-ray power spectrum has been 
estimated in the past (Tagliaferri et al 1991, Zhang 
ct al 1999) to have a slope of ~ —2, similar to the 
slope of the EUVE power spectrum of Mrk 421 
. Both objects are High frequency peak BL Lac 
objects, and their EUVE and X-ray emission are 
thought to be due to synchrotron radiation. The 
similarity in the power spectrum slopes could im- 
ply that in addition to the same emission mecha- 
nism, the same variability mechanism operates in 
both objects in the two bands. 4 

In general, red noise power spectra can result 
from a shot-noise model, i.e. from a superposi- 
tion of random "flare" -like events, each one with 
the same duration and shape, and with an ampli- 
tude randomly distributed around a mean value. 
A slope a = 2 is expected in the case of expo- 
nentially decaying shots, consistent with our best 
fit slope value. However, the shape of the flares 
in the EUVE light curves of Mrk 421 can by eye, 
be more complicated than that; they are asym- 
metric and they do not have the same shape all 
the times. For example, the first flare in the DS 
1995 lightcurve (Figure 6) shows a relatively slow 
rise (it reaches maximum after ~ 2 days), and 
faster decline (~ 1 day). Immediately after this 
the source rises again, reaching now the maximum 
after ~ 1 day. Even if it is not clear what happens 
after the maximum (if there is a second flare or 
not) this difference in rise times indicates that the 
flares are not always the same. Probably, there are 
flares which rise slow and decay fast and flares with 
the opposite behavior, or the complicated shapes 
can be the result of the superposition of consecu- 
tive flares. 



4 The normalized amplitude of the PKS 2155-304 X-ray 
power spectrum at 10 -5 Hz is > 1000 (Zhang et al. 1999), 
when for Mrk 421 the normalized amplitude of EUVE 
power spectrum at the same frequency, is 148 ± 33 (Ta- 
ble 2). 
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In general, knowledge of the power spectrum 
slope alone cannot provide us with much infor- 
mation on the variability mechanism that oper- 
ates in the source. One has to work the other 
way around; develop a model for the variability 
mechanism, predict the shape of flares, hence the 
power spectrum slope, and then compare it with 
the data. 

Even though more work is needed in order to 
investigate if any model can produce light curves 
with the correct power spectrum slope, amplitude 
and low frequency knee, it is interesting to notice 
the similarity between the shape of some of the 
flares in the EUVE light curves of Mrk 421 and 
those predicted by Celotti et al. (1991). The flare 
at the beginning of the DS 1995 light curve (Fig- 
ure 6) resemble by eye the flare expected in the 
case of a perturbation with growing thickness in 
an inhomogcneous jet (Figure 4b of Celotti et al. 
1991). A different behavior is present in the sec- 
ond flare in the same light curve; this time the flux 
of the source continues to be large after reaching 
the maximum and the amplitude of the event is 
smaller than in the previous case. This flare could 
correspond to the case in Figure 4a of Celotti et 
al. (1991), in which they used different parameters 
for the stationary jet structure. 

5.3. The low frequency turn over 

The presence of a low frequency knee, with a 
value that remains constant, within the error bars, 
in all the EUVE power spectra of Mrk 421 , indi- 
cates the presence of a "characteristic time scale" 
(t var ) in the system. This time scale determines 
the "memory" of the system, i.e. the variabil- 
ity mechanism operates in such a way so that, 
whatever the origin and nature of variations, they 
"last" typically for a period which is, on average, 
equal to t 

var ~ 1/(4 10 ^) Hz ~ 3 days. For ex- 
ample, according to the simple shot noise model 
mentioned above, the "memory" of the system is 
determined by the duration of the flares. If, on 
average, they last for t var days, a low frequency 
knee at — ^/t var is expected. As we mentioned 
in the beginning of this section, visual inspection 
of the Mrk 421 EUVE lightcurves implies that the 
flares in them last for ~ 3 days. 

The fact that the low frequency knee remains 
constant implies that t var is determined by a char- 
acteristic of the source that does not change with 



time. Such a permanent characteristic might be 
the size of the EUV emitting region. If that is the 
case, then an upper limit to the source size at rest 
frame, R, is given by R = t var x (5c) , where 5 is the 
Dopplcr factor. For t 

var ~ 3 days and a typical 
value of 5 ~ 10, the source size should be less than 
~ 8 x 10 16 cm. This value is similar to the value of 
4.8 x 10 16 cm and 4.7 x 10 16 cm that Chiaberge & 
Ghisellini (1999) and Mastichiadis & Kirk (1997) 
find from theoretical modeling of multifrequency 
energy spectra of Mrk 421 . 

Recently, Maraschi et al. (1999) proposed the 
e-folding times scales of 10.7xl0 4 s, 7.1 xlO 4 s, 
and 5.3xl0 4 satl,5 and 15 keV respectively for 
a large flare observed by BeppoSAX. These time 
scales were calculated from a simple arbitrary ex- 
ponential fitting to a single event in the BeppoSAX 
light curve. If they represent indeed a character- 
istic time scale of the system at these energies, 
then this time scale decreases with energy. The 
time scale derived from EUV data is 2.6 xlO 5 , i.e. 
larger than the time scales in the X-ray band. If 
the time scale is associated with the source size, 
then these numbers arc consistent with models in 
which the source size is expected to decrease with 
energy. A more detailed comparison has to wait 
for a proper statistical analysis of the X-ray data. 

An indication of a break at ~ 0.5 days is present 
also in the structure function analysis performed 
on BeppoSAX 2000 observation of MRK 421 
(Zhang Y. H. 2000 private communication). As- 
suming that the bulk of BeppoSAX emission is at 
-3-4 keV, the 6 times longer EUV (-0.1 keV) 
break could be explained with the synchrotron 
cooling of the same electron population emitting 
in the BeppoSAX range. 

5.4. Evidence of non-stationarity: DS 
1996 data 

As is evident from Figure 9 and Figure 12, the 
slope of DS 1996 power spectrum (open circles) is 
flatter than the others, showing the non station- 
ary nature of Mrk 421 emission. The appearance 
of the DS 1996 light curve (Section 2.2.4) agrees 
well with the power spectrum of that light curve. 
The flatter slope of the DS 1996 power spectrum 
indicates larger amplitude variations at high fre- 
quencies, i.e. lots of short time variations, consis- 
tent with the presence of a number of short lived 
flares, and the absence of any long "trend" in the 
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light curve itself. Changes in the power spectrum 
indicates that some of the variability mechanism 
characteristics change with time. PKS 2155-304 
also shows evidence for non-stationarity in its X- 
ray light curves (Zhang et al 1999). Mrk 421 has 
shown evidence of "non-stationarity" in the past: 
Fossati et al 2000, in a simultaneous TeV and X- 
ray observation, found that the hard X-rays were 
delayed with respect to the soft X-rays, opposite 
to the behaviour previously observed (Takahashi 
et al 1996). 

In the context of the shot-noise model, in DS 
1996 power spectrum: (a) the different slope in- 
dicates a difference in the flare shape; (b) the 
roughly similar amplitude (Figure 9) suggests that 
probably both the number of flares per unit time 
and their amplitude has not changed drastically 
during DS 1996 observation; (c) the knee at the 
same frequency, within the errors, of the other 
power spectra indicates that the duration of DS 
1996 flares is not changed. 

To the eye, though the flares in the DS 1996 light 
curve appear to be relatively short-lived. We spec- 
ulate that the DS 1996 flares consist of two parts: 
one with a fast decay/rise evolution, and another 
one with a much slower evolution, not seen di- 
rectly in the light curve, due to the blending of 
consecutive flares. 

6. Conclusions 

The extremely large EUVE database on Mrk 421 
(<~ 1100 ks) offers a so far unique opportunity to 
study the variability properties of a BL Lacertae 
object at high energies. We analyzed Mrk 421 
light curves by means of the power spectrum and 
structure function techniques; in a power spec- 
trum with slope —2.14 ± 0.28 we find indication, 
confirmed by the structure function analysis, of 
a break at ~ 3 days. This is the first time that 
a break in the power spectrum of a BL Lacertae 
object is found. We also find evidence of non- 
stationarity for Mrk 421 EUV emission. 

The few existing models on BL Lacertae objects 
variability (Celotti et al. 1991, Chiaberge & Ghis- 
cllini 1999, Georganopoulos & Marscher 1998), de- 
veloped for multifrequency analysis, try to explain 
the amplitude and delays observed among the var- 
ious energy bands. The single energy approach, 
as the one proposed in this paper, is at least as 



promising as the multifrequency approach, but no 
predictions on the power spectrum and on the 
structure function shapes (the most widely used 
analysis methods to quantify the variability of an 
object) have never been derived from any of the 
current models. 

It is beyond the scope of this paper to investigate 
in detail whether the model flares that have been 
computed so far can reproduce the results from 
the analysis of the EUVE light curves of Mrk 421 
and we reserve this for a future paper. However, 
such a quantitative work, will be very interesting 
since it may be able to distinguish between differ- 
ent models that have been proposed to explain the 
high frequency emission of blazars. 
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Fig. 1. — Comparison of the effective areas of the 
Deep Survey (solid line) and scanner A (dashed 
line) instruments on board the EUVE satellite. 

Fig. 2. — Simulated count rates in the DS (solid 
line) and Scanner A (dashed line) instruments as- 
suming an absorbed power law model with spec- 
tral index a = 2 The ratio of the count rate at the 
peak wavelength for the two instruments was used 
as rescaling factor to convert the Scanner A count 
rate into a DS count rate. 

Fig. 3. — Combined lightcurve of all EUVE obser- 
vations from 1994 to 1997 binned over one average 
EUVE orbit (~ 5544 s). 

Fig. 4. — April 2-12 1994 Deep Survey lightcurve 
binned over one average EUVE orbit (~ 5544 s). 
This light curve was also presented in Fruscione 
et al. (1996). 

Fig. 5. — February 4-7 1995 Scanner A light curve 
(open circles) binned over one average EUVE or- 
bit (~ 5544 s). The filled circles represent the 
count rate scaled to the DS instrument (see Sec- 
tion 2.2.2). 

Fig. 6. — Same as Figure 4 for the April 25 - May 
13 1995 DS data. The data from MJD 49832 to 
MJD 49843 were presented in Kartje et al. (1997) 
and were part of a multiwavelcnght campaign dur- 
ing wich a TeV flare was detected (Buckley et al. 
1996). 

Fig. 7.— Same as Figure 4 for the April 17-30 DS 
data (on the left) and the May 10-11 data (on the 
right). A large TeV flare (Zweerink et al. 1997) 
occured on 7 May 1996. 

Fig. 8. — Same as Figure 5 for the February 7-11 
1997 Scanner A data. 

Fig. 9. — Normalized power spectrum (binning 
factor m = 10, see Section 3.1) obtained from each 
lightcurve binned over 5544s (1 average EUVE or- 
bit). 1995 DS/S observation was split into two 
parts to avoid the gap. The open symbols repre- 
sent DS/S observations (triangles for 1994, squares 
for the first part of 1995, stars for the second part 
of 1995 and circles for 1996) while the filled sym- 
bols represent scanner observations (squares for 
1995 and circles for 1997). 



Fig. 10. — Same as Figure 9 for DS 1996 observa- 
tion with different binnings. Opens circles, as in 
Figure 9, represent a binning value of m = 10, and 
filled circles m = 20 (see Section 3.1) 

Fig. 11. — Average normalized power spectrum 
of Mrk 421 obtained from all the 5544s binned 
lightcurves presented in Figure 3 (see section 
3.1). Different symbols represent different bin- 
nings: m = 20 (filled circles) for the low fre- 
quency part and m = 10 (open circles) for the 
high frequency part representing the the Poisson 
noise level. The solid line shows the best power 
law model fit with slope -2.14. 

Fig. 12. — Power spectra of the DS 1995 (a) and 
DS 1996 (b) light curves binned over 500 s (open 
circles). The long dashed line in both plots indi- 
cate the expected Poisson noise power level. Solid 
line shows the best fit of the model defined by 
equation 5 in Section 3.2, with the slope fixed to 
-2.14. The dotted line in (a) represents the simu- 
lated power spectrum using the same best fit val- 
ues but moving the break to 10~ 7 Hz; while the 
dotted line in (b) shows the best fit of the same 
model but with the slope left as a free variable (see 
text). 

Fig. 13.— Same as Fig 12 but for DS 1994, scan- 
ner A 1995 and scanner A 1997 respectively. 

Fig. 14.— Structure Function of the total EUVE 
light curve of Mrk 421 . (Figure 3) Filled circles 
correspond to the SF estimated with a lag bin of 
size 1 day and the filled squares to the SF esti- 
mated with a lag bin of size 1 month. Open sym- 
bols (circles and squares) correspond to the mean 
SF of 1000 simulated light curves with a power 
spectrum of slope 2.14 and a frequency turnover 
at 4 x 10~ 6 Hz. Dotted triangles correspond to 
the mean SF of 1000 simulated light curves with a 
same slope power spectrum but with a frequency 
turnover at 10 -7 Hz. In both cases, the error- 
bars correspond to the standard deviation of the 
simulated SFs around their mean (see Section 4 
details). 
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Table 1 

EUVE Observations of Mrk 421 



Year 


Start 


End 




Instrument 3, 


-i' b 

J- cxp 

(ks) 


Count rate 6 
(10~ 2 c s- 1 ) 


Overall Variability 
factor 


1994 


2 Apr 


12 Apr d 




DS 


280 


18 ±0.85 


2 


1995 


4 Feb 


7 Feb 




ScaA 


68 


32 ±2.9 
45±4.1 d 


1.5 


1995 


25 Apr 


13 May 




DS 


355 


29 ±0.96 


3 


1995 


25 Apr 


28 Apr 


(flare) 6 






37±2.3 




1995 


29 Apr 


6 May 


(decay) 6 






29 ± 1.5 




1995 


7 May 


13 May 


(quiescence) 






22 ± 1.0 




1996 


17 Apr 


30 Apr 




DS 


299 


30 ± 1.08 


1.4 


1996 


10 May 


11 May 




DS 


3.6 


40 ± 1.1 




1996 


7 Feb 


11 Feb 




ScaA 


108 


19 ± 1.6 
27 ±2.2 


2 



a "DS" indicates that the source was observed in the Deep Survey photometer. "ScaA" indicates that the source 
was observed in one of the three EUVE photometers (scanner A) during a pointing within the Right Angle Program. 

b Total exposure time calculated eliminating (i) all SAA passages (ii) all satellite daytime data (hi) all times during 
which the detector was turned off (iv) all times affected by possible earth blockage. It does include corrections for 
telescope vignetting, deadtime and limited telemetry allocation (primbshing) . 

c Average count rate in the Lexan/B filter (« 60 — 180 A). For the scanner observations, the second line gives the 
count rate normalized to the one in the DS assuming an average spectral shape, to take into account the difference 
in the effective area between the two instruments. 

d This observation was presented in Fruscione et al. (1996) 
6 This observation was presented in Kartje et al.( 1997) 
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Table 2 

Power spectra "power law" model fits for the DS orbit-binned light curves of MRK 421 



Year 


Instrument 


a 


A 


No. of points 


x 2 


1994 


DS 


4.02 ± 1.42 


13 ± 24 a 


3 


0.38 


1995 (part 1) 


DS 


4.14 ±0.86 


22 ± 15 a 


4 


0.35 


1995 (part 2) 


DS 


2.64 ± 1.10 


15±8 a 


3 


0.17 


1996 


DS 


1.45 ±0.42 


13±6 a 


5 


3.85 


All 


DS+scaA 


2.14 ±0.28 


148 ±33 b 


10 


5.6 



a in units of (counts/sec) 2 /Hz 
b in units of Hz -1 
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Table 3 

Power spectra "power Law plus low frequency turn over" model fit for the DS and 

SCANNER A, 500 SEC BINNED LIGHT CURVES OF MRK 421 



Year Instrument Slope2a fb DoF % 2 

(Hz) 



1994 


DS 


2.14 (fixed) 


8.4tf 7 


x 10" 


-6 


43 


67.8 


1995 


ScaA 


2.14 (fixed) 


5.0tf 3 


x 10" 


-5 


14 


28.7 


1995 


DS 


2.14 (fixed) 


2-6+iig 


x 10" 


-6 


77 


116.9 


1996 


DS 


2.14 (fixed) 


7 7+4.3 
'•'-2.6 


x 10" 


-5 


53 


140.4 






(1.1 ±0.1) 


(1-6±1?b 


x 10" 


- 6 ) 


(52) 


(65.5) 


1997 


ScaA 


2.14 (fixed) 


™ + -l 7 s 


x 10" 


-5 


17 


16.7 
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